Process and apparatus for producing colored chemical coatings

ABSTRACT

A computer-aided color formulation process is provided wherein a color standard is matched by a blend of pigments from a library of pigments. Reflectance values of the standard color at a plurality of reflectance points are supplied to a computer programmed to select pigments from the library and to approximate concentrations of the selected pigments to give a match to the reflectance data of the standard. The computer is also programmed to generate the pigment composition from a group of those selected to match the color of the standard, and to determine the metameric index of the match. The color of a blend of library pigments prepared in response to the concentration values generated by the computer is measured by a colorimeter under daylight to insure that the difference in color coordinates between the blend and the standard is within tolerance. If the color coordinate difference is not within tolerance, the process provides for determining the amount of selected pigment addition necessary to shade the blend to provide a color coordinate difference within tolerance. The blend is then measured by the colorimeter under incandescent light to insure that the metameric index is within tolerance.

Unit l s tent [72] Inventor William H. McCarty lllunterdon, NJ. 859,4

Sept. 19, 1969 Aug. 24, 1971 Mobile Oil Corporation [2 l App]. No. [22] Filed [45] Patented [73] Assignee [54] PROCESS AND APPARATUS ll fllli PRODUCING COLORED Cl'lEMlCAlL COATIINGS 25 Claims, 3 Drawing Fig. [52] 235/150 Glloll 7/00 [50] 356/175; 235/150 [56] References Cited Inspect ColorOt Standard ll And Measure Reflectance Run Computer Program To Obtain Pigment Selection Program To Decide 0n Pigments Examine Pigment Selection And Decide On Pigments To Use Primary Examiner-Malcolm A. Morrison Assistant Examiner-Edward J. Wise Attorneys-Oswald G. Hayes, Andrew L. Gaboriault and James F. Powers, Jr.

ABSTRACT: A computer-aided color fonnulation process is provided wherein a color standard is matched by a blend of pigments from a library of pigments. Reflectance values of the standard color at a plurality of reflectance points are supplied to a computer programmed to select pigments from the library and to approximate concentrations of the selected pigments to give a match to the reflectance data of the standard. The computer is also programmed to generate the pigment composition from a group of those selected to match the color of the standard, and to determine the metameric index of the match. The color of a blend of library pigments prepared in response to the concentration values generated by the computer is measured by a colorimeter under daylight to insure that the difference in color coordinates between the blend and the standard is within tolerance. If the color coordinate difference is not within tolerance, the process provides for determining the amount of selected pigment addition necessary to shade the blend to provide a color coordinate difference within tolerance. The blend is then measured by the calorimeter under incandescent light to insure that the metameric index is within tolerance.

M ls Run Computer Program MI /7 20 To Obtain Concentrations k? ForColor Match Under Make Up Run mp ogram To ObtalnAmounts Of Daylight/ind Under ln m, p T Add candescent Light (Ml) lgment o /5 Measure N A 22 23 Difference In I O g fil s l-r ece In 51, i a? 6 2 Color Under 24 Wme Up y g Y Incandescent FOtmUlO FT /9 e5 m PATENTEU M1824 19.: 3,601 58E) SHEET 2 0F 3 Standard CO Spectrophometer Encoder ".33 Panel X 1 Computer 34 t Flow Rate l Controllers For m 35 I Pigments i t 'r ScreW Feeders 36 For Pigments 5 Color A Corrections 37 S olvents, gg? Resms,Etc.

Undler Computer Control w To Color Color Product Correctlon Testmg Storage Vessel vessel INVENTOR.

lPli tUClESS AND APPARATU lFOlR PRODUCING COlUDlltlElDt CHEMICAL COATINGS BACKGROUND OF THE llNVENTlON 1. Field of the invention This invention relates to process and apparatus for color formulation. More particularly, it relates to computer-aided process and apparatus for formulating colored chemical coatings.

2. Description Of the Prior Art Color is the psychological response of the observer to the sum of the visible radiant energy, reflected from an object. The factors that determine color are the natures of the light source, the eye of the observer, and the reflecting surface.

A light source can be characterized by measuring its intensity as a function of wavelength in the visible region of the spectra. Standard light sources are defined in Color in Business, Science, and Industry by Judd and Wyszeclri, Wiley, New York 1967), pp. 108. These are:

Source A incandescent light Source B noon sunlight Source C average daylight.

lf colors which match under one light source do not match under a different light source, the match is termed metameric." In general, a color match should have as little metamerism as possible.

The most difficult aspect of color is the nature of the eye of the observer. Since it is rare that two observers are identical, a standard or average must be defined. The standard observer used today was recommended in 1931 by the International Commission of illumination, and is discussed at page 111 of Color in Business, Science, and Industry noted above.

The standard observer is defined in terms of the amount of red, green, and blue primary lights (called the tristimulus values) required to produce a color match to light of a given Wavelength and unit intensity. The tristimulus values of the standard observer can be combined with the light intensity of the standard sources. This result is shown in Color in Business, Science, and Industry at pages 132-433. These results are called the tristimulus values of the standard sources A, B, and C,

The appearance of the surface of an opaque object is characterized by its reflectance spectrum which is a function of the wavelength of the incident light and is measured (generally at 10 millimicron intervals) with a spectrophotometer.

The tristimulus values for a surface illuminated by a standard light source can be calculated as the product of the reflectance and the tristimulus values of the standard source at each wavelength.

The color is the sum of the tristimulus values of the surface over the visible wavelengths. Since there are three tristimulus values at each wavelength, the result of the summation is three numbers, called the tristimulus coordinates of the color. A plot of all conceivable tristimulus coordinates generates a space in three dimensions called color space.

The quantitative measure of color coordinates is performed by a colorimeter which automatically measures the reflectance spectrum and performs the necessary multiplications and summations.

The problem in color matching is to generate an experimental surface that has the same tristimulus coordinates as the standard surface. This can be done by generating a surface which has the same reflectance spectrum as the standard. Such a surface would be a spectral match" to the standard. A spectral match is sufficient to gua a nonmctameric color match. That is, the surface generated will match the standard under any light source. A spectral match is not necessary to obtain a color match under a single light source. Surfaces having quite different reflectance spectra can be color matched. However, the matches will be metameric.

It has long been a problem in the chemical coatings art to match a standard color with a blend of pigments. For example, a coatings manufacturer may receive a sample standard from an automobile manufacturer with a solicitation for a bid on a predetermined amount of chemical coating having a color matching that of the standard. The chemical coatings manufacturer is then faced with the problem of matching the standard with a combination of pigments.

A trial and error method for color matching provides for inspecting the color of the standard and then guessing the pigments in a library and relative concentrations necessary to provide a color matching blend. A blend of pigments is then prepared based on the pigment and concentration guesses. The blend is sprayed on a sheet of metal to provide a panel and the panel is compared with the standard in daylight. If the visual inspection indicates that the standard and the panel do not match, the colorist guesses at the amount of pigments to add to the blend and the new blend is applied to another panel. The colorist again visually compares the standard and the panel and this procedure is repeated until the colorist visually observes a match. The colorist then observes the standard and the matching panel under incandescent light. lf there is no match under incandescent light, the colorist repeats the foregoing steps until a match is obtained under incandescent light. Thereupon, the colorist prepares a formula indicating the pigments and concentrations thereof necessary to provide a blend of pigments having a color which matches that of the standard.

Thus, in the trial and error method, colorists are required to exercise human judgement based on visual observation and experience.

It is also known to use a spectrophotometer for measuring reflectance values of the sample, and guess the pigments and pigment concentrations of the sample based on an inspection of the measured reflectance values. Further, it is known to use a colorimeter to determine the differences in color between the pigment blend and the sample under light sources A and C. However, it is still necessary to exercise human judgment.

SUMMARY OF THE lNv'ENTlON in accordance with the present invention there is provided a process of selecting pigments for a pigment library having a standard pigment to color match a colored surface coating. The process comprises operating an automatic processing system to perform the steps of generating signals representative of the reflectance spectrum of the surface coating, and generating signals representative of absorption constants and scattering constants of N library pigments N l. in response to these signals, the process generates signals representative of hl-l concentration ratios between nl nonstandard library pigments and the standard library pigment, in accordance R the corrected reflectance of the surface coating at a In accordance with another aspect of the present invention the concentration ratio C lC values are refined by nonlinear least square solution procedures based on the total differential of the dependent variable.

and generates signals representative of the weight fractions of the nonstandard library pigments considered from the library in accordance with l l/ l)( l)v where C, weight fraction of the standard pigment,

N the number of pigments considered from the library,

and

C /C the concentration ratio of the i"' pigment to the standard pigment.

In accordance with yet another aspect of the present invention, a signal representative of a predetermined weight fraction is generated. The weight fraction signals of the nonstandard library pigments are compared with the predetermined weight fraction signal, and the compared weight fraction signals having a value less than the predetermined weight fraction signal are eliminated from further processing. The process is then repeated until each of the remaining nonstandard pigments has a weight fraction at least as great as the predetermined weight fraction.

A further aspect of the invention provides a system for selecting pigments from a pigment library having a standard pigment to color match a colored surface coating. The system comprises means for generating signals representative of the reflectance spectrum of the surface coating, and means for generating signals representative of absorption constants (K) and scattering constants (S) of N library pigments, N 1. The system also includes means responsive to these signals for generating signals representative of N-l concentration ratios between N-l nonstandard library pigments and the standard library pigment, in accordance with:

where R the corrected reflectance of the surface coating at a given wavelength, K S the absorption and scatter constants of the standard library pigment at a given wavelength, K,-, S the absorption and scatter constants of the i'" pigment in the library at a given wavelength, N the number of pigments to be considered from the library, and C,/C the concentration ratio of the i"' pigment to the standard pigment. Thus, the present invention provides a process and a system for minimizing the need for exercising human judgement in selecting pigments to color match a standard color.

BRIEF DESCRIPTION OF THE DRAWINGS FIG. 1 is a block diagram of a computer-aided color formulation process,

FIG. 2 is a block diagram of a continuous process for blending a chemical coating, and

FIG. 3 is a block diagram of a general-purpose digital computer.

DESCRIPTION OF SPECIFIC EMBODIMENTS An Iterative Procedure for Solving Systems of Nonlinear Equations Since linear equation theory is the basis of nonlinear solution procedures used in this invention, an outline of the applicable linear theory is presented before proceeding to the technique used to solve nonlinear equations.

Assume that a theoretical relationship exists between a dependent variable y, and a number of independent variables X (i=1, 2...N) as given by equation (1).

The unknown coefficients, k,, are calculated by use of p different measurements of the dependent variable y, made at p different settings of the independent variables. The experimental measurement is designated by g}; (i=1, 2...p), and there results a system of p equations in N unknowns.

The system of equations can be written in matrix notation as given by equation (2),

Y=I 1 where 'y= a column vector whose elements are the predicted values ofy,(i=l ...p)

I a column vector whose elements are the unknown k,

(i=I N) [X] the matrix of independent variable settings, whose elements are X the subscripts designating the row and column respectively.

If p=N, the desired values of k, can be obtained by inverting X. Thus,

[X the transpose of [X] That is, the solution of equation (4) minimizes'S, where S is given by equation (5).

where The equations (2) and (4) are true for linear systems. However, they can be used as a starting point in the iterative solution of some nonlinear systems.

For example, suppose that a nonlinear relationship exists between y and x of equation (6), and that equation (6) can be rearranged to give equation (7 That is, equation (6) can be linearized so that the left-hand side contains only y and the right side is linear in the unknown k If the same set of experimental measurements were available as before, the set of equations (7) could be written in matrix notation as givep by (8).

{X'} matrix whose clcmcnts arcjly t Equation (8) resembles cquation (2.), except that the isle merits of the X matrix depend on the predicted value of the dependent variable. A. formal solution can be obtained if the elements of the X matrix are: evaluated using the experimental values of 31. Thus a first estimate oi lrcan be obtained by use of equation (9) o jli 'l l 'l WlJlJ l The value oil I: obtained may be us 10 predict the values of y by means of equation (6). in this way, a i, st estimate of the 1': and the predicted y, can be tounrl.

Successively better approximations to the desired least square solutions for in are obtained by utilizing standard nonlinear least squarc solution technique host: techniques are based on the total rlififercntial of the dependent variable, given by equation l).

D ral y 1 P 1- u ctr- "'1 J1 tani J t t where the (O l /r) In) are obtained by differentiating {t5 The set of equations i0) can be written in matriur notation as given by equation l i .Uy-" E Pjrilt EF-G t h 17?) Where y, isgbtained from (6) using the current values Gift That is, Dy is a vector that. points in a direction towards the experimental values of y. a

Substituting the P matrix for the X matrix and U32 iorfty -i in equation (9 leads to solution n'lt. The elements of die are the changes that must be added to the current estimate of into obtain an improved estimate. Thus, the new values of loan: given by equation l 3).

neyv old dki The new lnarc used to predict new valucs oty and the entire procedure is repeated to obtain a still better estimate of in. if successive iterations leacl to ptO' t cly smallct' chart In, the process is said to convert and may be terminated when the kmo longer change appreciably.

Since the solution obtained minimitcd the sum of the squared deviation between the experimental and predicted y, the procedure can be thought of as minimizing the distance squared (and thus the distancc) behvecn thc coordinates of the experimental and predicted y in a p dimensional space.

The liubelltalvlunk analysis described in tinfoiin Business, Science, and Industry by Judd and lll vszcclci, Wile New "torlt (1967 relates the light a ring and absorp propert .i of pigment dispersions to the observed rcl'lcctar c of a surface. lf the pigment dispcrr on d substrate, then the reflect the overall reflectance of tilt:

by lfiubelka and lvlunlr trc complete hiding.

Before proceeding to the escription oi the several assumptions oi the These assumptions are:

a. The surface reflectance lium having the same refractive index as the p .mcnt ve clc.

b. The pigment particles scatter li ht diff 'ly. Surface reflectance measurements are ma rly violating the first assumption. The eff ol the a rurtace interface is accounted for through the use of correction equation as given by l4).

m/(fi r two cases, or should he pointed out.

where R he reflectance measured in air R the equivalent reflectance measurement made in a medium of the same refractive index as the pigment vehicle.

The constants oz, ,6, 'y were determined at a color laboratory of the assignee of the invention as follows:

[1! 0.005 /Ei==0.57ti 'y=0.4 The second assumption of the theory is valid for spherically shaped pigments, but is invalid for tlaltelike particles, for example, flalted aluminum.

a. incomplete Hiding of the: Substrate The reflectance from a surface at incomplete hiding of the substance is given by equation where R the reflectance at incomplete hiding. R the reflectance at complete hiding.

Rg' the reflectance of the substrate Z=XZ1C;S;(1 R' )/R it the weight of pigment per unit area.

C 5 the weight fraction concentration and scatter constant of the i" pigment in the mixture.

[if the number of pigments in the mixture.

Equation 15) may be solved for N Argos.

as given by 16).

b. Complete Hiding oi the Substrate At complete hiding of the substrate, the relation between reflectance, pigment concentration, and pigment constants is given by equation 17),

where it,- the absorption constant of the i'" pigment and all other terms retain their former significance.

Equation (17) may be solved for R, to give l8) R'=1+awa +2a (1 where x EK i=1 5: N

The value of& depends on the relative values of the pigment constants, rather than their absolute values as demonstrated by dividing the numerator and denominator of 6 by 6], with the result shown in equation (20).

1 I) N s 9) =2 1 1 A IN A QIQ 2. Color Coordinates and Metamerism The tristimulus coordinates of a surface under a given light source are defined by equations (2l (22) and (23):

700 X K HXRAX 700 Y K H -RA)\ 700 Z= K H,RA)\

x 400 (23) where H,, H, H, the tristimulus values of a standard light source such as given in Color in Business, Science, and Industry, supra,

A the wavelength AA the wavelength increment at which the reflectance R is measured.

The coordinates X, Y, Z are transformed into the cube root coordinates L, a and b as given by equations (24 (25) and (26) The physical significance of the coordinates are as follows: L lightness a redness-greenness b yellowness-blueness The reflectance spectrum must be available to calculate the color coordinates of a surface. The reflectance may be measured on an existing panel, or the reflectance, R, may be calculated by use of the Kubelka-Munk equations (15) or (18) and corrected to give R by use of equation 14).

The difference in color between two panels is represented by the difference in color coordinates. If the desired or standard color is designated by subscript s, and the color of the current batchis given by subscript b then the color differences is provided by equations (30), (3 l and (32).

The color is matched when AL=Aa=Ab=0. Matching is done for the standard source C. When a match has been obtained, the values of AL, Aa, and Ab under source A are calculated to-define the index of metamerism, Ml, as given by equation (33).

The metameric index is an indication of the light source dependency of the color match. In general, it is desirable to find a match that has the smallest metameric index.

A high metameric index is due to a poor spectral match. Poor spectral matches occur when the standard and batch colors are made from different pigments.

3. Vectors and Shading The color coordinates (L, a, b) of a panel are functions of the concentrations of the pigments. That is L=L o) (34 a=a (G) (35) b=b (6 6) where Similar expressions can be written for the total differentials of a and b.

if a column vector dL were defined as having the elements dL, da, and db, then the change in color coordinates with respect to a change in concentration coordinates can be written in matrix notation as dL=[V]dC' (38) where 92 2A P 6C ac. OC 0a 0L 56, 'T 2 .25 60 DO The matrix V is the matrix of the transformation from differences in concentration space to difference in color space. In practice, the individual partial derivatives are called vectors." This terminology should not be confused with the notion of vectors in matrix algebra.

Equation (38) is strictly true only for infinitesimal changes. However, in practice, it is desired to match colors which are close to one another, but observably different.

For this reason the AL, Aa, Ab values representing the difference in color space between the standard and batch panels can be substituted into equation (38) to give the more practical equation (39). AL: [l'.]AC

Equation (39) may be used to estimate changes in pigment concentration to effect a color match. This is called SHAD- ING." However, the elements of the V matrix (Vectors) must be known.

Vectors may be obtained in two ways. They may be estimated experimentally by varying the concentration of the pig ments one at a time, and subsequently measuring the new color coordinates by use of a colorimeter. The ratio of the observed change in color coordinates to the change in concentration give the desired vectors.

Altemately, vectors may be calculated by use of the Kubelka-Munlt equation and the definition of the color coordinates. lFor example,Z)L/ DC, can be calculated by means of the chain rule as shown by equation (40).

acfac or OR DR Ta T.

j indicates one of 31 wave lengths at which R is measured by a spectrophotometer, and

iindicates the i pigment, i=1, 2...N.

All of the terms in equations (41 (42) and (43) can be calculated, given the set of K and S values for each of the pig' ments in the batch.

Assuming that the vectors have bee obtained and that the difference in color coordinates between standard and batch has been calculated, then equation (39) can be solved for the AC vector necessary to effect a color match. However, there are only three color coordinates. Thus, only three pigment composition changes can be calculated. As discussed hereinabove under Complete Hiding of the Substrate, if N pigments are present, only N-l of these are independent. These considerations lead to the following rules for the use of equation (39), depending on the value of N.

I. If N=3, only two pigments can be varied. in this case, no exact color match can be obtained. However, the least squire solution will minimize the sum of the squared deviation between the standard and batch colors. That is, the distance in color space between the standard and batch colors may be minimized. If the matrix, V, contains only the vectors corresponding to the two pigments to be varied, then the solution is given by equation (44). A

It is immaterial which two pigments are selected.

2. If N=4, then three pigment compositions can be varied and AC is given by equation (45).

ill

In this case the Ml is invariant. The matrix V contains only the vectors corresponding to the pigments selected. The solution is independent of which three pigments are chosen.

3. If N 4, then there exists (Nl )l/3!BN-4 )!combinations of three pigment composition changes that may be calculated by means of equation (45 Hence, there are an equivalent number of color matches. Not all of these matches will have the same Metameric Index however. The best color match will have the smallest Ml.

In any event, the new concentration vector is 7|Ell 0d where C,,,,,is the current value of C, and AC contaigs the changes calculated via equations (44) or (45); i.e., AC will either have two elements or three elements depending on the number of pigment&

The elements AC are of cgurse added to the correspondi n g elements in C to obtain C The concentration v ector C may be normalized by dividing each element in C by the sum of the elements in C Of course, the relationship between pigment composition and color coordinates is nonlinear, while equations (44) and (45) are based on linear equation theory as discussed hereinabove under An Iterative Procedure for Solving Systems of Nonlinear Equations. Thus, after one application of equation (44) or (45), the new color coordinates for the batch should be calculated. It may be found that the new batch color coordinates are still not the same as the standard coordinates. Thus, equation (41) or (42) must be used repeatedly, each time using the color coordinates and vectors based on the previously predicated pigment concentrations.

If the standard and batch colors are sufficiently close to start with, this iterative procedure will converge; i.e., AL, Aa, and Ab will approach zero or constant values with continued iteration. However, if the initial colors are too widely separated, the method r r iay diverge. In this case, the change in concentration vector, AC, is too large. To avoid divergence, AC, should be multiplied by a fraction less than one before calculating the new concentration.

The following Table I lists the nomenclature of the equations used in the present specifications and appended claims.

TABLE I NOMENCLATURE a color coordinate in L, a, b system b color coordinate in L, a, b system B color coordinate in G, R, 8 system C fractional composition D the total differential of the quantity F indicates a functional relationship G color coordinate in G, R, B system H,, H H, tristimulus values for the standard light sources K absorption constant; normalization factor in x, y. z coordinate system k an unknown L color coordinate in L, a. b system Ml metameric index N the number of pigments R color coordinate in G, R, B system;

reflectance data point or reflectance spectra R reflectance spectra in medium of same refractive index as the vehicle of the coating at complete hiding of the substrate R" same as R except incomplete hiding of substrate S scattering constant [V] matrix of partial derivatives relating concentration space to color space X the red tristimulus coordinate of a standard light source; the pigment weight per unit surface area [X] a matrix of independent variable settings y a dependent variable Y the green tristimulus coordinate ofa standard light source lnbit' l unttnucd Z the blue tr'istimulus coordinate of a standard light source Greek Symbols A wavelength in millimicrons a difference =a partial derivative C/IO 6 defined by equation l9) Superscripts a vector 0 a measured value a modified value of the variable superscript as defined in the text T transpose of the quantity -l the inverse of the quantity subscripts PROGRAM OPERATION FIG. 1 shows a block diagram representation of a computeraided color formulation process in accordance with the present invention. Following the start 10 of the process, a colorist inspects the color of the standard which is to be matched by a blend of pigments from his library of pigments, and measures the reflectance of the standard color at a plurality of reflectance points. The reflectance values are fed to a computer programmed to select pigments from the colorists library and to approximate concentrations of the selected pigments to give best match to the reflectance data of the standard color 12.

Signals representative of the selected pigments and concentrations thereof are automatically fed to a program for choosing a subset of the selected pigments 14. Four pigments are eferably initially chosen by the program 14. Alternatively the computer provides a readout to the colorist representative of the selected pigments and the concentrations thereof. The colorist then decides on a number of the selected pigments 13, preferably initially four, for use in the following program 15.

Generally, the program 14 and the colorist decision 13 use those pigments having the greatest approximate concentrations generated by the selection program 12.

In response to signals representative of the measured reflectances and to signals representative of the chosen pigments in program 14 or to the pigments chosen by the colorist 13, a formula program 15 is entered to obtain concentrations of the chosen pigments necessary for a color match under daylight to generate signals representative of the metameric index (Ml) of the match. The MI is automatically compared with a predetermined Ml tolerance, or the Ml is read out to the colorist for comparison with the predetermined Ml tolerance 16.

If the metameric index (Ml) is unsatisfactory, the program 14 is automatically reentered and another subset of the pigments selected by the selection program 12 are used, if such program selected a number of pigments greater than the number ofpigments initially chosen by the program 14.

Alternatively, if the Ml is unsatisfactory, the colorist again examines the pigment selection from the program 12 and chooses another subset of the pigments selected by the selection program 12 for input to the formula program 15. The loop including the program 14 and the program 15 or the loop including the colorist decision 13 and the program 15 are repeated until a satisfactory Ml match is obtained 16. If the MI is still unsatisfactory after repeating the loop including the program l4 and the program 15 until all subsets of the pigments selected by the selection program have been inputted to the program, it is necessary for the colorist to then exercise human judgement based on his experience 13 to decide on the set of pigments with which to enter the program 15.

When a satisfactory Ml match is obtained 16, a blend of pig ments is prepared 17 in response to output signals from the program 15 which are representative of the pigments and concentrations thereof for providing a satisfactory Ml match.

The difference in color between the standard and the blend 17 is measured by a colorimeter under daylight 18, If the color coordinate differences between the standard and the blend are not within tolerance 19, a program 20 is entered to shade or vary the concentrations of the pigments in the blend 17. This loop is repeated until an acceptable color coordinate differences tolerance is met 19.

When an acceptable color coordinate differences tolerance is obtained 19, the standard and the blend 17 is measured by a colorimeter under incandescent light to obtain a measure of the metameric index (Ml). If the MI is not within tolerance, the process returns to program 14 or to the colorist decision step 13, and the process is repeated until a measured acceptable Ml tolerance is obtained 22 between the standard and the pigment blend. After an acceptable Ml tolerance is obtained 22, the colorist prepares a formula 23 representative of the library pigments and concentrations thereof necessary to provide a color match with the standard, and the process thereafter is ended 24.

It has been experienced that when the formula program 15 computes a color match having a satisfactory metameric index, the verification steps 18, 21 usually confirm the computed value of metamerism within acceptable limits of error. Thus, the verification steps 18, 21 may be eliminated, and production runs carried out based on the pigment concentration values generated by formula program 15. The production blend may then be sampled and measured by a colorimeter under daylight as described with reference to block 18. Errors in pigment concentration may then be corrected by theshade program 20.

A suitable instrument for providing reflectance values in step 11 is a spectrophotometer model 505 manufactured by Bausch and Lomb, Rochester, N.Y. The reflectance values may be automatically punched on paper tape by a Datex encoder manufactured by Datex Corporation in the form of integer numbers between 001 and 999, representing fractional reflectance values between 0.001 and 0.999. The instrument may be operated at a gain of 10, so that reflectances between 00l and 0l0 may be amplified by a factor of 10. Before running a panel having a standard color thereon on the spectrophotometer, the number of decimal places (either three or four) is punched manually on the tape to identify the gain level. The data representative of reflectance values is than punched automatically by the spectrophotometer at predetermined millimicron intervals between 400 and 700 millimicrons. A suitable interval is 10 millimicrons to provide a total of 31 reflectance points. The reflectance values at the 3l points together with the gain level identification constitute a reflectance spectrum. The reflectance spectrum may be read directly into a computer through a Teletype such as manufactured by Teletype Corporation, Skokie, lll. A suitable colorimeter for use in measuring color differences 18, 21 is the Colonnaster manufactured by Manufacturers Engineering and Equipment Corporation, Warrington, Pa.

The following is a discussion of a specific embodiment described with reference to programs which appear in attached appendices which form part of the present application. The programs are computer listings written in BASIC. A discussion of the BASlC language appears in BASIC Langauge: Reference Manual (71 1224A), General Electric Company and the trustees of Dartmouth College, 1968. The pro grams were prepared for a time-sharing service provided by' General Electric Company. The service and system commands are described in Mark II Time-Shan'ng Service Command System: Reference Manual (711223), General Electric Company, 1967. A computer suitable for the service is the GE. 635. Each of the foregoing manuals 711223 and 71 1224A are incorporated herein by reference.

The first three programs SCATTER, KS1 and KS2 are used to generate absolute pigment constants for all library pigments in a given product line.

As discussed hereinabove with reference to equation (16), the standard pigment is applied to a plurality of experimental panels each at different levels of incomplete hiding of the substrate. A panel is also prepared with the standard white applied at complete hiding of the substrate. The reflectance spectrum of each of the panels is then obtained for use in program SCAT'llER.

An additional plurality of experimental panels are prepared each having a chemical coating of complete hiding of the substrate and containing a different concentration of a plurality of pigments from the library including the standard pigment. The reflectance spectrum of each of the panels is then obtained for use in programs KS] and KS2 PROGRAM SCATTER The SCATTER program which appears in Appendix 11 is used to calculate absolute scattering and absorption constants for the standard pigment in a product line which is defined by the vehicle, e.g., alkyl or acrylic. in this program a white pigment, W812, was chosen as the standard pigment. The program uses equations (14), l6), and 17).

For one pigment, equations (16) and (17) reduce to (47) and (48) The experimental data are entered into a SDATA file as shown in Appendix 1. The reflectance spectrum at complete hiding is entered first, followed by the reflectance spectrum of the substrate. The data at incomplete hiding is entered next, with the weight of pigment per sq. ft. of surface entered immediately after the number of decimal places in each spectrum. The system commands and data entered by the user is underlined.

Execution of program scatter provides the result shown in Appendix III. The user must specify the file name to which K and S values are to be written and the number of panels at incomplete hiding that have been entered for each substrate entered. The program computes the scatter values, at each wavelength, for each panel in the substrate set using equations (14) and (47). At the completion of the calculation, the average S value for all panels in the substrate set is calculated. The program then repeats the calculation for the next substrate set, if there is another. If not, as in the example considered, the program stops when the SDATA file is completely read in.

The program will operate with any number of panels per substrate set, and any number of substrate sets. The program output is the value of S at each wavelength, averaged over all reflectance spectra at incomplete hiding in the SDATA file.

Using the grand average values of S, the program computes the corresponding K values by use of equation (48). After printing these values at the terminal, the program writes the K and S values to the file designated by the user, Appendix IV shows a listing of the W812 file which was written by the program.

PROGRAM xsl Program KS1 shown in Appendix V provides a more precise measure of the K values of the standard pigment, and a measure of how well a given set of data can be fitted by the Kubelka-Munlt theory.

The program will calculate the K values of the standard pigment, and the 1C and S values for all other pigments in mixtures containing up to six pigments. Up to 40 reflectance spectra may be used in the calculation.

The data are entered into the DATA file as shown in Appendix VI. The concentration data is entered first, for each panel, with the concentration of the standard pigment given first. The corresponding reflectance spectra are then entered in the same order as the concentration data.

Execution of the program is demonstrated is demonstrated in appendix VII. The user is asked for the file name of the pigmerits, the number of pigments and the number of panels. in the example considered, there are four pigments and 21 panels. The pigment codes of the four pigments used are indicated in Appendix VII.

The calculation proceeds by first obtaining an initial guess of the values of the constants. This is done with the aid of equation (49), which is a rearrangement of equation (17), in a form that is linear in the unknown constants.

N E as,

where subscript 1 indicated the standard pigment.

The S values of the standard pigment, 8,, the concentration, and the reference data are used to evaluate the left-hand side, and the coefficients of the unknown K, and S,- in equation (49), at each wavelength.

The least square solution procedure outlined hereinabove under Iterative Procedure for Solving Systems of Nonlinear Equations is followed to obtain the least square estimate of the constants. The predicted R are then calculated via use of equations l8) and (14).

The nonlinear least squares procedure outlined hereinabove under iterative procedure is then followed to obtain the least square values of the constants.

The wavelength, the sum of squared deviations and the mean square error are printed out at the terminal. A is then printed, to which the user is instructed to reply with a l or a 2 depending on the amount of printout desired. Illustrations of the printout for both responses is given in Appendix VII.

The calculation proceeds, wavelength by wavelength, until wavelength 700 is reached, at which time the values of K and S for each pigment are written to files, and execution ter minates. The calculated K and S values are shown in Appendix Vlll.

PROGRAM KS2 where subscript 1 indicates the standard pigment.

Equation (50) is the same as (49), with the first term in the K summation moved from the right-hand side of the equation to the left-hand side. The concentration and reflectance data are placed in the DATA file as in Appendix Vl.

Program execution is illustrated in Appendix X. The user is asked for the number of pigments, the number of panels and the pigment codes of the pigments used in the mixture. The program uses the K and S values of the standard pigment, the concentration. and the reflectance data to evaluate the lefthand side of equation (50) and the coefficients of the unknown K. and 8,. After obtaining the initial estimates of the constants, the program proceeds to calculate the nonlinear least square solutions for the unknown K, and S, according to the procedure described hereinabove under the iterative procedure.

Following the calculation at each wavelength, the program prints the wavelength, the sum of squared deviation, and the mean square error between experimental and predicted reflectance. The user may request a more detailed printout at each wavelength by typing a 2 in response to the printed after each wavelength calculation.

After completing the calculation at wavelength 700 the program writes the calculated K and S values for each pigment to files bearing the same name as the pigment codes supplied by the user. The pigment constants calculated by KS2 are given in Appendix Xl. These results are based on the values of K and S given in Appendix 1V, rather than those shown under W812 in Appendix VIII.

The pigment mixtures used in the experimental panels need not contain all the pigments. It is recommended that the masstones (single pigment dispersions) of each pigment, except the standard, be included in the experimental panels. The remaining panels may consist of a number of single pigment dispersions let down in the standard, as well as mixtures of all the pigments.

The advantage of being able to evaluate a number of pigment constants simultaneously is to provide pigment constants that fit reflectance data ranging from single pigment dispersions, to mixtures of eachpigment in the standard, to mixtures of a number of pigments, and to determine the accuracy of the fit to all these different types of dispersions.

Thus, the absolute K and S values for the library pigments may be determined by executing program SCATTER to generate K and S values of the standard library pigment. Program KS1 is entered with the S values of the standard as determined by program SCATTER and with concentration and measured reflectance data to generate improved K values of the standard pigment and K and S values for up to five nonstandard library pigments. Program KS2 is then entered with the improved K values of the standard pigment from program KS1, the S values of the standard pigment from program SCATTER, the known pigment concentration values, and the measured reflectance data to generate K and S values for the remaining nonstandard library pigments.

Alternatively, program KS2 may be entered directly with the K and S values of the standard library pigment from program SCA'ITER, the known concentration values, and the measured reflectance values to generate the K and S values for the nonstandard library pigments.

lf absolute K and S values are not desired, relative values may be obtained as described hereinabove under Complete Hiding afthe Substrate.

PROGRAM SELECT Assuming that the pigments K and S values have been obtained, the problem of selecting pigments from the library to match a given color of unknown composition remains to be solved.

The problem is one of qualitative and quantitative analysis. It is necessary to know which pigments to use and how much of each pigment should be used. The available data consists of the 31 reflectance values of the unknown. This data represents the combined effect of all the pigments (generally 3 to 5) actually present in the panel.

In accordance with an important aspect of the present invention, the measured reflectance values of the unknown or standard panel, to be matched by a blend from the library, are used to select those pigments from the library that best fit these data. This aspect comprises using a linearized form of the Kubelka-Munk equation (17) to arrive at a first estimate of the concentration of pigments in the library necessary to give a best match to the reflectance spectrum of the standard panel.

Examination of equation (50) indicates that if one divides through by the concentration of the standard white, C the resulting equation is linear in the concentration ratio C,/C as shown by equation (51 5 N [351' (K.+s.) .]g S.))Z-Q1+R s. so.

where R the corrected reflectance of the standard at a given wavelength.

K S the absorption and scatter constants of one pigment chosen as the standard (a white pigment is best), at a given wavelength.

K S, the absorption and scatter constants of the i'" pigment as the library, at a given wavelength.

N the number of pigments to be considered from our library.

C,/C the concentration ratio of the i" pigment to the standard pigment in the library.

The reflectance data of the standard panel measured in block 11 of FIG. 1 and the library pigment constants are used to evaluate the left-hand side and the coefficients of C,/C in equation (51) to obtain a first estimate of C,/C, via the technique described hereinabove under Iterative Procedure for Solving Systems of Nonlinear Equations.

The values of C,/C may be further refined by use of the standard nonlinear least square solution procedure of equations (l0) through (13).

Since there are 31 reflectance values per unknown, and there are N-l independent concentration ratios, up to 32 pigments may be included in the calculation.

A scan of the estimated concentrations of each pigment indicates the relative amount of each pigment. By dropping out all pigments present in amounts less than some small positive value, a number of pigments can be eliminated. The process can be repeated, using only the remaining pigments, until all pigments remaining are present in amounts greater than the amount chosen.

Program SELECT shown in Appendix XIII is based on equation (51) and is designed to perform the scanning and elimination of pigments described above.

The panel code, the product line, the number of decimal places and the panel gloss are punched on paper tape, followed by the 3l reflectance values for the panel. Panel gloss will be discussed hereinafter with reference to program FOR- MULA. The tapes are then entered in the RDATA file as shown in Appendix Xll. All of the K and S values are stored in the select program along with the corresponding pigment codes.

The execution of program SELECT is shown in Appendix XlV. The user is asked for the reflectance set, the tolerance and the number of pigments. The meaning of these terms are:

a. Reflectance Set the position of the desired panel in RDATA file.

b. Tolerance the weight fraction composition below which pigments are to be thrown out of the mixture.

c. Number of pigments the number of pigments to use initially in fitting the reflectance data (the number of pigments are taken in the order in which they are stored in the program).

3,601,589 ll? lid Program SELECT reads the pertinent data from the lrihltll unllnucd RDATA file, prints out the code, the line. the number of decimal places. and the gloss. The program then proceeds to g back cglnn ngo e re. calculate the initial estimate of the pigment concentration F 320430 Give line feed lamina, and prim ratios and the corresponding weight fraction composition of J H hcadingirtformation Fss.z7.wt. each pigment according to equations (52) and (53 and beint the P product line number. the number of decimal places and the panel gloss respectively. 1 I 340 Move the decimal place. Wl places to 61* 1 C C1 the left. for all reflectance data in T the spectrum. 342-348 Correct reflectance data for gloss as C a 6.(C pereq;uation(54).

x a 350-37!) Store the integer num ers 1 through N5 in vector Q. 380 Set iteration counter N to I. Program SELECT scans the weight fraction concentration, 04 set matrices to proper size depending and throws out any pigment present in amounts less than a value 075 6430-450 Give line feed and print out the tolerance specified. The number of pigments under consideration is printed out, along with the pigment codes thrown out. This procedure is repeated as many times as necessary, to ar- 4 7 rive at a set of pigments all of whose composition exceed the tolerance, at which time the answer is printed out and the user is asked to specify a new tolerance. If the new tolerance is the same as the original, the program skips to the beginning and asks for the next reflectance set, the tolerance and the number of pigments. 580.630

Table ll provides an explanation of Program SELECT as presented in Appendix Kill].

TABLE ll e60 [00-110 Designates disc file name where the reflectance spectra of the standards are stored and designates that individual data points are separated by blanks. l20-l50 Serve to designate the maximum size of the matrices used in the calculation. [-180 Define the values ol'the constants used in the conversion of R to R in equation (I4).

I90 Read the number 5 which is the 4 5 800 number of sets of pigment constants stored in the program(this causes the value 24 to be read from line 1090 320 200-220 Read the pigment code names stored in lines 1100 through H30 into the 5 vector P5.

230 Read all the pigment constants into matrix file read the data beginning on line 1140.

240 Causes a line feed at the telctypewriter.

250 Print the statement "Ref. Set.

Tolerance. No. of Pig." This alerts 5 5 the user that certain information 10204050 will be required as follows: Causes a break in program execution. 10604080 Program will wait for the values of l T. and N5 to be input by the user.

N] the position ofthe 6O reflectance spectrum in the RDATA file. for

which a pigment number of pigments. 5 currently under consideration.

Correct the reflectance data calculate R from R according to equation (14); calculate the left-hand side of equation (5 l and store in vector B; calculate the coefficients of the unknown pigment concentration ratios. C,/C.. according to equation (51) and store in matrix A.

Calculate the unknown CJC, ratios according to the least square solution procedure as described in iteration procedure. and summarized by equation (9).

Set a counter L to 5 Set a counter M to I.

Go to subroutine beginning on line 860 and calculate the estimated weight fraction of each pigment. store results in vector W and return to line 670.

Scan the elements in vector W. and

except the first drop out all pigments present in amounts less than T. except for N801. Test N801 for negative value ifiteration counter is I; this block ofcodc results in a redefinition of the vector 0. such that 0 contains the subscnipt numbers ol'the pigments retained by the test. The counter L contains the number of pigments whose subscripts appear in Q.

lncrcment the iteration counter by l.

A remark statement summarizing the result of the operations defined by lines 670-790.

Test to see ifthe number of pigments after the scan (namely L) is the same as the number before the scan. If so. print out the pigment codes and estimated concentrations according to lines 960 through 1010. if not set N5 to L and go to line 390 and repeat the entire calculation for the new number of pigments.

Set 11 T, and ask user for new tolerance.

If new tolerance. T. is same as previous (stored. in .11) then go to line 240 and repeat entire procedure for another case. Otherwise. using the new tolerance. return to line 670 and res-can the weight fraction concentrations.

selection is desired.

T the tolerance. or

weight Traction concentration. below The results of six trial runs with program SELECT using which pjgmcmamw panels of known composition are given in Table III. These be dropped. results clearly show that the correct pigmentation is selected.

.\'5-= the num er of pigments to consider in but that several spurious pigments may be retained.

The pigment codes used in Table III refer to commercially the selection (this Wm .hnaluc of available pigments, the specific identification'of which are not N5 read initially). necessary for understanding the invention. To practice the in- 17040 Rcad heading "Mmammnd vention, the pigments in the library must be nonflakelike partireflectance data from the RDA'JA file. NI times. The result is to place the reflectance data of the N1" cles which sufficiently scatter light diffusely to satisfy the second assumption discussed hereinaboveunder Kubelkapanel in matriit L. MunkAnalysis.

TABLE Ill.-(().\ll.-\RIFON F SELE TED \llTll ACTUAL IIUMENTATION ute l Actual Pitrinmt liizmcnt Panel .\'0. L'Odl' ('unv. code Cour.

R 11 20,1 W81! 0. $054 :81? 0. $506 X5 1 0141' \S01 0. 015" L507 0521 LQOT 0. 0617 (7216 0567 G210 0. 0118 L-SU'J U, 0302 L50 0. 0104 YSta 0. 0144 R 3 (2 13 U. 6011 W 12 0. 655 X801 0. 0054 \501 0v 0040 YSlb 0. 3704 1 515 0. 333 G215 0. 00776 U216 0. 0065 YSIT 0. 00737 R 33115) X801 0. 00311 5:801 0.0028 YSlli 0. 2277 t 818 0. 232 [515 0. 7106 Y815 0. 7651 R 44 ('20) W812 0. 818 W81? 0. 8542 X501 0. 0165 \'S01 0. 0154 G216 (I. 033! G216 0. 0203 L800 0. 0371 L806 0. 020 Y417 0. 0045 Y41T 0. 0811 R 55 (21) W812 0. ['6 W81? 0. 072 L807 0. 00746 L807 0. 0058 Y807 0. 0127 1'80? 0. 0104 120 0. 010'.) Y120 0. 0126 R 66 (24) W812 0.151 W812 0. 102 N501 0. 0026 \801 0, 0027 Y815 0. 111 Y815 0. 120 YSIS 0. 60 Y818 0. 776 YSOT O. 04

' Number of pigments considered in the selection. I 2 W is white, N is a black. L is a blue. G is a green, Y a yellow, and R i a red.

Extensive tests have shown that one or more subsets of the pigments chosen by the SELECT program will result in a satisfactory color match in about 70 to 80 percent of the cases examined.

Thus program SELECT greatly reduces the prior art need for human judgment to choose the initial pigments.

PROGRAM FORMULA the FORMULA program shown in Appendix XV is designed primarily to compute the pigment composition necessary to obtain a color match, and the metamerism (Ml) of the match as described hereinabove with reference to block 15 of FIG. 1. The program has a number of other functions which will be discussed below.

Before discussing the program, a current listing of a typical pigment file is illustrated in Appendix XV. The first 62 values e the K and S values of the pigment. The 63rd value is the pigment cost in $/lb., the 64th value is the pigment density in lb./gal., and the 64th thru 67th value are the weight fraction composition of the pigment in the three tint paste qualities that are currently part of the color matching system.

The input data for FORMULA are entered in the RDATA file as shown in Appendix XII. The data are; the panel code, the product line, the number of decimal places, and the panel gloss, followed by the 31 reflectance values of the standard panel as measured in block 11 of FIG. 1.

All of these data have been discussed previously, with the exception of panel gloss. The color matching system was developed for a gloss level of 95 percent, as determined by a gloss meter. Therefore, any standard submitted must be corrected to the equivalent 95 percent gloss level. This is done thru the use ofequation (54).

R=R,,,(95-G) (0.028)/l00 (54) where R,,, the measured reflectance R the reflectance corrected to the 95 percent gloss level G the percent gloss Execution of FORMULA is shown in Appendix XVll. The

user is asked first for the starting reflectance set in the RDATA file. The program reads the RDATA file to the desired starting point, and print out the code, line, places and gloss which were read in from the file.

The user is asked for the number of pigments, and the pigment codes that should be used in the calculation. The program reads the pigment files, computes and prints out the (least square) spectral match composition, and the sum of the squared deviation between experimental and predicted reflectance. The program uses equation (51 )to obtain the initial estimate of the spectral match composition and then uses the nonlinear least square technique to arrive at the final composition. The compositions are in the same order as the pigment codes supplied.

The program proceeds to compute and print out the color match composition, the corresponding tint paste composition, the Metameric index, and the deviation between predicted and experimental color coordinates in the L, a, b system, under source C and source A. If the color match is perfect the differences in coordinates under source C will be zero. This condition is not desired in practice, since deviations less than 10 in color coordinates are not detectable.

The program proceeds to an option point, from which the user may exercise 9 options. These options are:

The use of the various options are illustrated in Appendix XVII. Options 1, 2, 3 and 8 are self-explanatory.

Option 4 of Appendix XVll illustrates the use of the program to reduce the difference in color coordinates under Source C to a lower value.

Option 5 of Appendix XVll gives an example of the print out for vectors. The row wise order of the vectors is 2a Ob ac, a0, so,

The vectors are used to shade the match in the laboratory after it is made up and the actual color coordinates are determined. Shading proceeds as described hereinabove under Vectors and Shading with the exception that changes in pigment composition must be positive. It can be shown, that for any pigment mixture at least one set of three vectors will give concentration changes that are all positive.

Option 6 in Appendix XVll shows the input and output information for hiding. The user is asked for the substrate number, the color tolerance, the pigment volume percent (PVC), the dry film thickness (DFT), and the weight fraction solids (WFS) in the wet paint. The significance of these input data are as follows:

Substrate Number designates the position, of the reflectance spectrum of the substrate which is to be hidden. These data are placed in a file called SUB.

Color Tolerance The square root of the sum of the squared deviations in the color coordinates, L, a, b under Source C between the color at complete and incomplete hiding that is desired.

PVC

the percent by volume of pigment in the dry paint.

DFT

the dry film thickness which is to be used to hide the substrate to the tolerance specified. WFS

ments are kept out for tinting, and to uncontrolled variations from batch to batch.

Pigment additions are computed as described hereinabove under l/eclors and Shading with the exception that the vectors the weight fraction of pigment plus resin in the wet paint. remain constant for each iteration, and the pigment concen- In specifying these input items, the user cannot specify the tration changes must always be positive because pigments cantolerance, and PVC, and DFT. Only two of these quantities not be removed from a batch. As discussed previously, it can are independent. In order to show which of these three items be shown that there exists a set of positive composition are to be computed, the user enters a zero for the unknown changes for at least one combination ofthree pigments. quantity. An example of the operation of the program is illus- The following Table IV provides a summary of the functions trated in Option 6 in Appendix XVII. and limitations of the foregoing programs.

TABLE IV Program name Function Limitations SFA'ITE IL... Calculatc absolute absorption and scatter constants for the standard pigment None. KS1 Calculate improved values 01 the absorption constant for the standard pigment using Up to mixtures containing up to 6 rcllcctancc spectra taken from mixtures of pigments with the standard. pigments. KS2 Calculate absorption and scatter constants for pigments mixed with the standard pig- D0.

mcnt. SELECT v Select pigmcnts from library for use in formulating a color of unknown composition Up to 32 pigments may be present in the library (Generally one more than the number of reflectance data points.) FORMULA. Given a set of pigments calculate: Up to 7 pigments may be used.

1. The pigment composition necessary to match a color. 2. The metamcrism of the mate 3. The pigment volume percent, or the dry film thickness, necessary to hide a surface of known reflectance, and the pigment cost per 100 sq. ft. of surface. 4. Vectors or use in shading laboratory formulations. I 5. An initial load formula and the cost of the pigment per gallon of wet paint. \ECTO R Given an existing color formulation, calculate: Up to 6 pigments may be used.

1. The color coordinates.

2. Correct thc formulation to the color coordinates of the standard, if necessary. 3. Calculate tint paste vectors for use in shading plant batches.

Option 7 of Appendix XVII illustrates the print out of the reflectance difference between the batch and standard colors for the spectral match compositions.

Option 9 of Appendix XVII illustrates the operation of the program for writing a summary of the calculation to a file. The user may be matching colors submitted from a variety of locations. By naming files for the submitting locations, the user may write a summary to the appropriate file for later recall and transmission to the submitter. Appendix XVIII shows the output for the two different formats used.

. PROGRAM VECTOR Program VECTOR used to shade in block 20 of FIG. l is shown in Appendix XX. The input data to program VECTOR are entered in the VDATA file as shown in Appendix XIX. The code of the color, the line, and the number of pigments are entered, followed by the pigment codes and concentration for each pigment or tint paste present in the formula.

Execution of the program is shown in Appendix XXI. The user is asked for the starting Data set. The corresponding code, line, number of pigments, and the pigment composition are read from the file and printed at the terminal.

The program computes the L, a, b coordinates for the color, under source C and asks the user to check the values against the L, a, b of the standard. If the L, a, b of the formula agree with the standard, the user responds with a yes" and l percent tint paste vectors are printed out. The columnwise order of the vectors is the same as the order of the pigments, and the rowwise order is bL/5C,,Oa/HC, ,bb/i'ic The vectors represent the changes in L, a, b coordinates per pound of tint paste, for an infinitesimal addition to a 100 pound batch.

The program output for the case where the L, a, b coordinated do not agree with the standard is illustrated in Appendix XXI. In this case the user is asked to supply the L, a, b of the standard. The difference in coordinates is then computed by the program, and the output compositions are shaded to the coordinates of the standard. When this process is completed, the vectors are printed out, followed by a revised formula which has the same weight of pigment as the original.

This feature of the program is used primarily to bring existing product color formulations into line with the instrumental and computer color program.

The vectors obtained apply at the L, a, and b of the standard. Actual plant batches differ from the standard by significant amounts, due in part to the fact that a portion of the pig- FIG. 2 schematically shows a continuous process for blending a chemical coating. With reference to FIG. 2, at the start 30 of the process a standard color panel 31 is placed in a spectrophotometer 32, e.g., the Bausch and Lomb-505, for determining the reflectance values of the standard. The reflectance values are automatically punched on paper tape by an encoder, e.g., the Datex Encoder. The data on the tape is fed to a computer 34, e.g., the GB. 635. The: connection between the encoder 33 and the computer 34 may comprise a Teletype (not shown) and a telephone line (not shown) as described in Mark Il-Time-Sharing Service: Reference Manual (711223), General Electric Company, 1967.

The computer 34 executes program SELECT to select library pigments and concentrations thereof to match the reflectance spectrum 32. The computer 34 then performs the routine described with reference to block 14 of FIG. 1 to decide on the initial set of library pigments to be considered. Program FORMULA is then executed, and if necessary the routine of block 14 of FIG. 1 is repeated until an acceptable metameric index is generated.

The computer 34 then generates output signals for adjusting library pigment flow rate controllers 35. The flow rate controllers 35 control screw feeders of tanks of library pigments 36 such that the library pigments flow at rates sufficient to provide concentrations for a color match as generated by program FORMULA.

The library pigments flow to a mixing vessel 37 along with solvents, resins and other components as indicated by the arrow entering at the right side of the mixing vessel 37.

The output stream of the mixing vessel 37 is fed t o a testing vessel 38 where the blend is sampled for testing by a colorimeter 39, e.g., the Colormaster, to determine whether the color coordinate difference under daylight between the standard and the blend is within tolerance. If the difference is unsatisfactory, the computer performs a routine including program VECTOR to generate shading output signals shown by the dashed line. The shading signals control adjusts a duplicate set'of flow rate controllers as indicated by the dashed line in 35. These controllers control a duplicate set of library pigment screw feeders as indicated by the dashed line in 36. The duplicate set provides a color correction flow of library pigments to a color correction vessel 40 to shade the blend flowing thereto from'the color testing vessel 34. The output from the color correction vessel 34 is fed to product storage.

It is also contemplated that a batch process can be carrieo out in a manner similar to that shown in FIG. 2 with the flow rate signals being replaced by absolute values to pass predetermined amounts of library pigments to a mixing vessel. This vessel is then sequentially sampled for testing by the colorimeter to determine color coordinate difference in daylight and to determine the metameric index.

A block diagram of the general-purpose digital computer General Electric time-sharing system is illustrated in FIG. 3. The primary computer units are identified as input 41, central processing unit (C.P.U.) 42, storage 43 and output 44.

The input 41 is entered by a Teletype, e.g., Teletype Corporation Model 33ASR, and a telephone line 45. In accordance with the specific embodiments, the Teletype 47 is fed a tape output from an encoder (not shown). ln-this manner, the teletype 47 receives and transmits data, including measured reflectances and the programs. The input 41 receives the data and programs and transmits them through the central processing unit (C.P.U.) 42 and to the storage 43. The central processing unit 42 processes the data in accordance with the stored programs and transmits a readout to the Teletype 47 via the output 44 and a telephone line 46.

"A endix 11 (continued) TABLE V.PROGRAM CHARACTERISTICS SUMJMARY FOR GENERAL ELECTRIC MARK II SERVICE fiaize Maximum number of- CPU time I lrogram Primary function of program acters Pigments Panels Compile Run SCATTER Calculate absolute scatter and absorption constants for standard pigment. 2,000 1 100 0. 16 1 K I Recalculate absorption constant ofstandard pigment and constants for other 4, 000 6 40 0.35 55 pigment in mixtures.

KS2 Calculate pigment constants....... 6 40 0. 35 42 SELEC'P Select pigments from pigment library. 32 1. 70 3. 8 FORMULA Calculate pigment composition for color m h under Sou: 7 1. 25 5. 5 VECTO R. Calculate vectors for use in shading plant batches 6 O. 74 1. 2

1 Typical time (CPU means central processing unit of computer). 2 Not applicable.

. Appendix I lncer1n| Data in the SDAI'A HI:

iiii FILE NAME 5 DA READY J as READY TAP READY 9m REID Y B52 READY Appendix II Program Scatter I00 FILES SDATMVSIZ Appendix II(continued) LET .G-B TB 200 Appendix III Executic'n of Prgmnm SCATIEII SCATTER 09:59 07/02/69 FILE NAME FER CBNSTARTST P0812 110. 01' PANELS I11 SUBSTRATE SET? 3 BF PIGHEMT/ SQUARE FEET: 0.00142 5 VALUES FER P411151. 1 ARE 19 0F PIGHERT/ SQUARE F087: 0.002703 5 VALUES FER P111121. 2 ARE 6 VALUES FER PANEL 3 ARE AVERAGE VALUES 01" S 111 SET AVERAGE VALUES 8F S 111 ALL SETS VALUES BF 11 ARE WRITING VALUES 0? 11 4110 5 1'0 P 11812 FILE USED 1 .23 UNITS Standard White Pimgnt File Wriiten g SCATTER Annendix V A nendix V (continued) Program KS 1 REM PRGGRAM CBMPUTES PIG. CBHSTANTS FZR PIG. MIXTURES 220 MAT READ 1 1 ,C(112,N5)

92 REM CBNTAINING UP I?! 6 P16. USING UP TB 40 PANELS. 221 REM MAT R CDNTAINS THE REF. DATA 94 REM CGNSTANTS ARE RELATIVE T8 THE SCATTER VALUES UJBTAINED 222 FUR 1:1 T8 N2 98 REM FREM THE FILE NAMED FIRST BY THE USER. 223 READ A1 .21

100 FILES DATMMBIZ 65 224 mm READ mmum) DELIMIT 1'1 225 MAT M=( 10Y'Z1 WM 1 10 DIM A(40. 1 1 ),B(40, 1 .C(40,6) .D(1 1 ,40),,E(1 1 .1 l 226 FER J:l 78 31 DIM "'(11 ,1 1 ),G(1 1 1 ,H(1 1 ,1 1 I131 1 1 ),14(1 .40) 227 LET R(I,J)=M(1 ,J)

DIM P1 1 40) ,R(40,31 .T(10.40),U( 11 ,1 ),X( l 40) 228 NEXT J 135 DIM X( 1 ,31 229 NEXT 1 PR1 ST 118. BF PIG. 110. BF PANELS"; 240 PRINT INPUT 115,112 70 250 PRINT "PIG. CBDEMGIVE STANDARD FIRST)";

REM V2 IS THE WAVE LENGTH CQURTER 260 MAT INPUT P$(115) LET V2 =1 265 PRINT REM N1 IS THE NUMBER 0F CWNSTARTS 270 FILE #2,PS( l LET 171:2115-1 275 MAT READ 1 2,1 1 ,31

LET C5: .576 280 MAT READ 1'2,S( l .31)

200 LET 04: .005 282 MAT K: 1 Of-SMK 210 LET C5: .4 285 MAT S:( IOT-SN S 215 REM C CBNTAINS THE 0011C. DATA 75 300 MAT B=ZER(N2,1

141 501116111 VI Appen'dbt VII (cont'th) Entcrlnfi 05m in the DATA File D 01.0 211.1: NAME fiL 0,1,0 READY 01.0 1112 114112-- 5 1 1 11 11240Y 5 5g 11s1 15105 01/01 /05 112401 110. 01' P10. 110. 01' 141121.01 4.21 A Y l1%m P10. 00020101112 2141104110 1 1110111 PVSILPNBOIEPYSIS'PYSH TYPE 1 011 2 111 112901102 10 1 0.041 0.0000 0.000 0.0521 (M919 0.0009 0J0 on my 10 1111: 1 IF 041.0. nzsuus 4112 1101 02011120 2:22 2122 2111; 2422 4 4444 44 44444 1 420 1 012942-5 1 200212-0 1 1 0.2550 0.0010 0.5020 0.2001 0.2024 0.0041 0.4121 0.5202 223 E12 3 l 0.514 0.0041 0.5021 0.0120 450 1 '514102-5 1 '211012-0 1 T 0.5054 0.0041 0.2052 0.4241 60 232N755 55 6 11.4105 0.0041 0.1025 0.4105 15 4 1 5 D was 0 0099 o 3175 0 557B 410 2.015502-5 2.241412-0 1 T D'IOM 0'00); 0'2562 0'625 450 1 .145052-5 1 .54201 2-0 1 1' 0'05" 0'00" 0'20 o'nsz 400 2.014002-5 2.221102-0 1 1 500 2.500002-5 l .910002-0 1 I 0.0200 0.0020 0.1400 0.021 '0 n 5 webs 7 l 0.0102 0.0051 0.051 1 0.055 5 0.999 0.00083 0 o 520 2.50011 2-5 1.022552-0 1 l o 9907 o 0093 o o 500 2.001502-5 2.001012-0 1 i 09127 0'02 0 u 20 540 2.021102-5 2.02001 2-0 1 1 {9222 0 0 550 5.5222-5 2.100202-0 1 1 o 5 0 0 500 I v 4.121052-5 5.001512-0 1 o 5695 5 m5 0 o 510 0.405142-5 2.042012-0 1 1 I 018904 0 0 h 550 2.204052-5 0.011212-0 1 i 5 12 12 15 15 14 14 14 10 10 10 15 20 41 01 115 200 501 404 500 521 555 550 551 500 540 550 551 0 20 20 20 20 2o 20 20 20 22 24 21 20 22 20 50 55 52 100 101 505 401 410 511 555 542 544 544 540 541 545 545 5 21 21 21 21 2 20 20 22 51 54 50 40 40 41 41 05 110 200 021 420 400 522 500 545 545 540 541 540 541 545 2 02 52 55 4 00 50 4 4s 52 52 55 51 15 121 210 512 500 40 401 404 405 402 402 402 450 450 5 50 50 00 50 00 41 45 40 5 00 00 0a 00 0s 14 04 150 255 502 402 441 450 405 402 401 401 401 400 400 400 0 5 50 50 50 01 05 00 12 10 02 110 155 102 200 212 214 210 200 200 504 520 055 042 545 544 545 040 500 502 500 054 555 5 14 14 14 15 10 11 10 20 20 55 40 54 55 51 50 50 10 104 110 101 105 101 100 101 195 190 102 101 100 101 0 10 10 20 20 22 24 21 5o 00 5o 15 05 100 102 102 101 121 122 200 0112 455 405 514 522 524 524 525 525 525 521 520 0 15 15 10 10 10 11 20 25 02 41 10 04 104 100 100 100 120 102 212 251 201 214 215 210 210 215 211 210 200 200 204 5 12 12 12 15 14 10 10 24 02 41 50 115 152 151 142 205 202 505 520 551 040 541 541 550 005 555 555 501 520 5 11 11 11 12 15 15 10 24 52 51 00 154 102 205 200 215 225 240 210 205 200 202 202 201 205 205 204 202 210 211 2 0 0 424 502 005 041 044 044 040 052 052 051 050 040 041 045 041 041 045 042 042 042 040 040 040 000 000 054 050 000 050 050 000 5 204 505 512 512 012 512 510 514 010 512 011 510 500 505 501 000 505 005 002 502 250 200 201 205 204 202 200 250 201 200 220 5 155 100 100 110 110 110 110 110 110 100 101 101 100 104 104 105 102 101 100 150 150 151 151 155 154 152 155 151 140 140 140 5 12 15 15 15 14 14 14 14 14 12 12 11 11 10 1o 00 00 00 01 0'1 00 00 05 04 04 05 05 02 02 00 00 4 510 510 510 505 501 200 205 200 224 210 211 215 211 205 202 202 255 252 240 244 245 240 254 250 225 220 224 224 215 215 211 4 T 101 150 155 151 140 140 144 141 1 40 155 151 120 121 1 15 1 1 4 112 111 101 105 104 100 05 00 05 00 50 50 05 00 54 4:50 50 40 51 50 41 410 45 42 44 45 41 40 50 50 55 55 50 55 54 54 54 20 21 21 25 20 24 25 22 22 1124 Y n 1 11%111 Appendix v11 :38 1. 2 07g1E-4 ?.g52E-6' 31 l 50 9 -4 5 112-5 1 010 1.120012-4 1 .52512-5 1 T 0 20 2.01 1512-4 1 .541022-5 1 T 01:15 FILE NAME-1&1 65 050 2.545012-4 1 .005202-5 1 g 25?. I 11 1 01141511 1111 (00110 0.

; 0.521 0.520101 110. 112 FIG. 142. 01 241121.51 4.21 0.555 0.550500 70 0.550 0.502400 P10. 0002510112 0141104110 11110111 111212, 111201 111015,,PY015 .42; 0.451220 .4 0.400242 1Y1 2 1 011 2 111 1120201122 10 1 0.505 0.054105 W22 1 11 041.0. 112001.10 4112 1401 02s11120 0.000 0.55416 0.5 11 0.005502 1111112 1.2110111 sum 02 00041120 142411 0004112 21111011 0.044 0.000022 400 1 .101212-5 5.5112-1 '1 1 75 0.100 0.200040 410 1.5-04542-5 1 .040052-0 1 1 0.522 0.520100 

1. A process of selecting pigments from a pigment library having a standard pigment to color match a surface coating, wherein said library pigments are nonflakelike pigments which scatter light diffusely, comprising the steps of: generating signals representative of the reflectance spectrum of the surface coating; and controlling an automatic processing system in response to a stored program to perform the steps comprising: a. generating signals representative of absorption constants (K) and scattering constants (S) of N library pigments, N>1, b. in response to said reflectance spectrum signals and said absorption and scattering constant signals, generating signals representative of N-1 concentration ratios between N-1 nonstandard library pigments and the standard library pigment in accordance with: where R'' the corrected reflectance of the surface coating at a given wavelength, K1,S1 the absorption and scatter constants of the standard library pigment at a given wavelength, Ki,Si the absorption and scatter constants of the ith pigment in the library, at a given wavelength, N the number of pigments to be considered from the library, and Ci/C1 the concentration ratio of the ith pigment to the standard pigment, c. and generating, in response to the N-1 concentration ratio signals, signals representative of library pigments selected for use in matching the color of the surface coating.
 2. The process of claim 1 wherein said stored program controls the automatic processing system to perform the additional step of: refining the concentration ratio signals obtained in step (b) by nonlinear least square solution procedures based on the total differential of the dependent variable prior to performing step (c).
 3. The process of claim 1 wherein the concentrations Ci and C1 are weight fractions and wherein step (c) comprises the steps of: d. generating a signal representative of the weight fraction of the standard library pigment in accordance with: and e. generating signals representative of the weight fractions of the nonstandard library pigments considered from the library in accordance with: Ci (Ci/C1) (CL), i 2...N, where Ci weight fraction of the ith library pigment, C1 weight fraction of the standard library pigment, N the number of pigments considered from the library, and Ci/C1 the concentration ratio of the ith pigment to the standard library pigment.
 4. The process of claim 3 wherein said stored program controls of the automatic processing system to perform the further steps of: f. generating a signal representative of a predetermined weight fraction value, g. comparing the nonstandard weight fraction sIgnals with the predetermined weight fraction signal, and h. eliminating from further processing the compared weight fraction signals having a value less than the predetermined weight fraction signal.
 5. The process of claim 4 wherein said stored program controls the automatic processing system to perform the further steps of: i. repeating steps (a), (b) and (d) through (h) after eliminating from further processing the absorption and scattering constants of the nonstandard pigments having a weight fraction value less than the predetermined weight fraction value as determined in step (g) to thereby reduce n by the number of nonstandard pigments eliminated, and j. repeating step (i) until each of the remaining nonstandard pigments has a weight fraction signal with a value at least as great as the predetermined weight fraction signal.
 6. The process of claim 1 wherein said reflectance spectrum signal generating step comprises generating signals representative of reflectance values of the surface coating at a plurality of wavelengths to provide the reflectance spectrum of the surface coating.
 7. The process of claim 6, wherein the automatic processing system comprises storage means, and wherein said reflectance spectrum signal generating step comprises: measuring the reflectance values of the surface coating by a spectrophotometer, and entering the reflectance values in the storage means to provide a source for generating said reflectance value signals.
 8. The process of claim 1 wherein said reflectance spectrum signal generating step comprises generating signals representative of reflectance values of the surface coating at a plurality of wavelengths between 400 and 700 millimicrons to provide the reflectance spectrum of the surface coating.
 9. The process of claim 8 wherein said reflectance wavelengths are at 10 millimicrons to provide 31 reflectance values.
 10. The process of claim 1 wherein the automatic processing system comprises storage means, and wherein the standard library pigment absorption constant and scattering constant signals generated in step (a) are initially generated by the steps comprising: d. measuring the reflectance at a plurality of wavelengths (1) of each of a plurality of panels each having the standard library pigment as a chemical coating at a different level of incomplete hiding, (2) of a panel having the standard library pigment as a chemical coating at a level of complete hiding, and (3) of the substrates; e. storing signals representative of the measured reflectances in said storage means; f. storing signals representative of the weight of the standard library pigment per unit area on each of the panels in said storage means; and g. controlling said automatic processing system in response to another stored program to perform the steps comprising: h. in response to the signals stored in steps (e) and (f), generating a signal representative of a scattering constant of the standard library pigment in accordance with: where S scattering constant, R'' the corrected reflectance at complete hiding at a given wavelength, R'''' the corrected reflectance at incomplete hiding at a given wavelength, Rg'' the reflectance of the substrate at a given wavelength, and X the weight of the standard pigment per unit area, i. repeating step (h) at each wavelength at which the reflectance is measured in step (d), j. averaging the generated scattering constant signals in each substrate set to generate signals representative of the scattering constant of the standard library pigment at each wavelength averaged over all reflectance spectra at incomplete hiding, k. in response to the signals generated in steps (e) and (j), generating signals representative of absorption constants of the standard library pigment in Accordance with: K S(1-R'')2/2R'' where K absorption constant of the standard pigment, S scattering constant of the standard pigment, and R'' the corrected reflectance at complete hiding at a given wavelength, and (l) entering the standard library pigment scattering constants signals generated in step (j), and the standard library pigment absorption constants signals generated in step (k) in said storage means.
 11. The process of claim 10 wherein the nonstandard library pigment absorption constant and scattering constant signals generated in step (a) are initially generated by the steps comprising: m. measuring the reflectance at a plurality of wavelengths of each of a plurality of panels each having a chemical coating of a different concentration of a plurality of nonstandard library pigments and the standard library pigments; n. entering signals representative of the reflectances measured in step (m) and of the concentration of the library pigments used in step (m) in said storage means; and o. controlling said automatic processing system in response to still another stored program to perform the steps comprising: p. in response to the scattering constant signals entered in said storage means in step (l) and the signals entered in said storage means in step (n), generating signals representative of estimates of the absorption and scattering constants of the nonstandard library pigments at each wavelength in accordance with: where R'' the corrected reflectance of the surface coating at a given wavelength, S1 the scattering constant of the standard library pigment at a given wavelength, Ki,Si the absorption and scatter constants of the ith pigment in the library, at a given wavelength, N the number of pigments to be considered from the library, Ci the concentration of the ith pigment, and C1 the concentration of the standard library pigment, q. in response to the estimated scattering and absorption constant signals, generating improved nonstandard pigment scattering and absorption constant signals at each wavelength by nonlinear least square solution procedures based on the total differential of the dependent variable, and r. entering the improved nonstandard pigment scattering and absorption constant signals generated in step (q) in said storage means.
 12. The process of claim 10 wherein the nonstandard library pigment absorption constant and scattering constant signals generated in step (a) are initially generated by the steps comprising: m. measuring the reflectance at a plurality of wavelengths of each of a plurality of panels each having a chemical coating of a different concentration of a plurality of nonstandard library pigments and the standard library pigment; n. entering signals representative of the reflectances measured in step (m) and of the concentrations of the library pigments used in step (m) in said storage means; and o. controlling said automatic processing system in response to still another stored program to perform the steps comprising: p. in response to the scattering and absorption constant signals entered in said storage means in step (l) and the signals entered in said storage means in step (n), generating signals representative of estimates of the absorption and scattering constants of the nonstandard library pigments at each wavelength in accordance with: where R'' the corrected reflectance of the surface coating at a given wavelength, K1,S1 the absorption and scatter constants of the standard library pigment at a given wavelength, Ki,Si the absorPtion and scatter constants of the ith pigment in the library, at a given wavelength, N the number of pigments to be considered from the library, Ci the concentration of the ith pigment, and C1 the concentration of the standard library pigment, q. in response to the estimated scattering and absorption constant signals, generating improved nonstandard pigment scattering and absorption constant signals at each wavelength by nonlinear least square solution procedures based on the total differential of the dependent variable, and r. entering the improved nonstandard pigment scattering and absorption constant signals generated in step (q) in said storage means.
 13. The process of claim 1 further comprising the steps of: d. selecting a subset of the library pigments selected in step (c); and e. controlling said automatic processing system in response to another program to perform the steps comprising: f. repeating step (a) for the selected subset of library pigments to generate signals representative of an initial estimate of a composition of the subset to spectrally match the surface coating, g. in response to the subset composition signals generated in step (f), generating improved subset composition signals by nonlinear least square solution procedures, and h. generating a signal representative of the metameric index of the improved subset composition signals.
 14. The process of claim 13 further comprising the steps of: i. comparing the metameric index signal to a reference to determine whether the metameric index is within tolerance, j. when the metameric index is outside tolerance, repeat steps (d) through (i) with different subsets of the library pigments selected in step (c) until the generated metameric index signal is within tolerance.
 15. The process of claim 14 further comprising the steps of: k. when the metameric index signal is within tolerance, preparing a chemical coating composition in response to the improved subset composition signals, and l. determining the color coordinate difference between the surface coating and the prepared chemical coating composition.
 16. The process of claim 15 further comprising the step of: m. shading the prepared chemical coating composition with additional amounts of the subset of pigments when the color coordinate difference measured in step (l) is out of tolerance.
 17. The process of claim 15 further comprising the step of: m. measuring the metameric index of the prepared chemical coating composition.
 18. The process of claim 17 further comprising the steps of: n. repeating steps (d) through (m) with a different subset of the library pigments until the measured metameric index of the prepared chemical coating is within tolerance.
 19. A system for selecting pigments from a pigment library having nonflakelike pigments which scatter light diffusely including a standard pigment to color match a colored surface coating comprising: storage means for storing signals representative of the reflectance spectrum of the surface coating, and for storing signals representative of absorption constants (K) and scattering constants (S) of N library pigments, N>1; and processing means, programmed to operate on said reflectance spectrum signals and on said absorption constant and scattering constant signals in said storage means, for generating signals representative of N-1 concentration ratios between N-1 nonstandard library pigments and the standard library pigment in accordance with: where R'' the corrected reflectance of the surface coating at a given wavelength, K1,S1 the absorption and scatter constants of the standard library pigment at a given wavelength, Ki,Si the absorption and scatter constants Of the ith pigment in the library, at a given wavelength, N the number of pigments to be considered from the library, and Ci/C1 the concentration ratio of the ith pigment to the standard pigment, said processing means being further programmed to operate on the N-1 concentration ratio signals for generating signals representative of library pigments selected for use in matching the color of the surface coating.
 20. The system of claim 19 wherein the concentrations Ci and C1 are weight fractions and wherein said processing means is programmed to operate on the N-1 concentration ratio signals by: generating a signal representative of the weight fraction of the standard library pigment in accordance with: and generating signals representative of the weight fractions of the nonstandard library pigments considered from the library in accordance with: Ci (Ci/C1) (CL), i 2...N, where Ci weight fraction of the ith library pigment, C1 weight fraction of the standard library pigment, N the number of pigments considered from the library, and Ci/C1 the concentration ratio of the ith library pigment to the standard library pigment.
 21. The system of claim 20 wherein said processing system is further programmed for generating a signal representative of a predetermined weight fraction value, for comparing the nonstandard weight fraction signals with the predetermined weight fraction signal, and for eliminating from further processing the compared weight fraction signals having a value less than the predetermined weight fraction signal.
 22. The system of claim 19 further comprising: means for measuring the reflectance at a plurality of wavelengths (1) of each of a plurality of panels each having a standard library pigment chemical coating at a different level of incomplete hiding, (2) of a panel having the standard chemical coating at a level of complete hiding, and (3) of the substrates; means for storing signals representative of the measured reflectances, and for storing signals representative of the weight of the standard library pigment per unit area on each of the panels; and said processing means being further programmed to: a. operate on said measured reflectance signals and on said weight signals for generating signals representative of scattering constants of the standard library pigment in accordance with: where S scattering constant, R'' the corrected reflectance at complete hiding at a given wavelength, R'''' the corrected reflectance at incomplete hiding at a given wavelength, Rg'' the reflectance of the substrate at a given wavelength, and X the weight of the standard pigment per unit area, b. average the generated scattering constant signals in each substrate set to thereby generate signals representative of the scattering constant at each wavelength averaged over all reflectance spectra at incomplete hiding, and c. operate on said measured reflectance signals and on said averaged scattering constant signals for generating signals representative of absorption constants of the standard library pigment in accordance with: K S(1-R'')2/2R'' where K absorption constant of the standard pigment, S scattering constant of the standard pigment, and R'' the corrected reflectance at complete hiding at a given wavelength.
 23. The system of claim 22 further including: means for storing signals representative of the rEflectance at a plurality of wavelengths of each of a plurality of panels each having a chemical coating of a different concentration of a plurality of nonstandard library pigments and the standard library pigment, and of the concentrations of the library pigments; and said processing means being further programmed to: d. operate on said signals representative of the reflectance at a plurality of wavelengths of each of a plurality of panels each having a chemical coating of a different concentration of a plurality of nonstandard library pigments and the standard library pigment, on said signals representative of the concentrations of the library pigments, and on said signals representative of averaged scattering constants for generating signals representative of estimates of the absorption and scattering constants of the nonstandard library pigments at each wavelength in accordance with: where R'' the corrected reflectance of the surface coating at a given wavelength, S1 the scattering constant of the standard library pigment at a given wavelength, Ki,Si the absorption and scatter constants of the ith pigment in the library, at a given wavelength, N the number of pigments to be considered from the library, Ci the concentration of the ith pigment, and C1 the concentration of the standard library pigment, and e. operate on the estimated scattering and absorption constant signals for generating improved nonstandard pigment scattering and absorption constant signals at each wavelength by nonlinear least square solution procedures, based on the total differential of the dependent variable.
 24. The system of claim 22 further including: means for storing signals representative of the reflectance at a plurality of wavelengths of each of a plurality of panels each having a chemical coating of a different concentration of a plurality of nonstandard library pigments and the standard library pigment, and of the concentrations of the library pigments; and said processing means being further programmed to: d. operate on said signals generated by (b) and (c), on said signals representative of the reflectance at a plurality of wavelengths of each of a plurality of panels each having a chemical coating of a different concentration of a plurality of nonstandard library pigments and the standard library pigment, and on signals representative of the concentrations of the library pigments for generating signals representative of estimates of the absorption and scattering constants of the nonstandard library pigments at each wavelength in accordance with: where R'' the corrected reflectance of the surface coating at a given wavelength, K1,S1 the absorption and scatter constants of the standard library pigment at a given wavelength, Ki,Si the absorption and scatter constants of the ith pigment in the library, at a given wavelength, N the number of pigments to be considered from the library, Ci the concentration of the ith pigment, and C1 the concentration of the standard library pigment, and e. operate on the estimated scattering and absorption constant signals for generating improved nonstandard pigment scattering and absorption constant signals at each wavelength by nonlinear least square solution procedures based on the total differential of the dependent variable.
 25. A process of selecting pigments from a pigment library to color match a colored surface coating, wherein said library consists of nonflakelike pigments which scatter light diffusely, comprising the steps of: generating signals representative of the reflectance spectrum of the sUrface coating; and controlling an automatic processing system in response to a stored program to perform the steps comprising: a. generating signals representative of absorption constants (K) and scattering constants (S) of N library pigments, N>1, b. in response to said reflectance spectrum signals and said absorption and scattering constant signals, generating signals representative of N-1 concentration ratios between N-1 nonstandard library pigments and the standard library pigment in accordance with a form of the following equation which is linear in the concentration ratio: where R'' the corrected reflectance of the surface coating at a given wavelength, Ki,Si the absorption and scatter constants of the ith pigment in the library at a given wavelength, N the number of pigments to be considered from the library, and Ci the concentration of the ith pigment, c. and generating, in response to the N-1 concentration ratio signals, signals representative of library pigments selected for use in matching the color of the surface coating. 